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Within the phenomenological approach a criterion for a choice of 
£5 the order parameter of the superfluid phases of 3 He in aerogel in a 

| vicinity of the transition temperature is derived. The order parame- 

Ch ' ter of bulk B-phase of 3 He meets this criterion and that of the bulk 

A-phase doesn't. A class of order parameters of Equal Spin Pairing 
(ESP) type meeting the derived criterion is discussed. Order parame- 
^vq | ters belonging to this class are proposed as candidates for the observed 

t> ■ A-like superfluid phase of liquid 3 He in aerogel. Effect of magnetic field 

, on the order parameters of this class is considered. 
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Among the other physical objects with unconventional Cooper pairing the 
*T3 • superfluid phases of 3 He are best understood. The pairing is referred as 

unconventional when simultaneously with the gauge symmetry some other 
q ! symmetries of a normal (nonsuperconducting or nonsuperfluid) phase are 

broken. Liquid 3 He, filling empty space in a high porosity aerogel can be 
superfluid [1]. The aerogel in that case plays a role of impurities. This 
system provides a good possibility to investigate effect of impurities on the 
superfluidity of 3 He and more generally on unconventional Cooper pairing. 
This knowledge can be applied to metallic unconventional superconductors 
UPt 3 , UBei 3 , Sr 2 Ru 2 4 , UGe 2 etc., where different types of impurities are 
inevitably present. Aerogel is different from usually considered point-like 
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impurities. It can be modelled as a rigid frame consisting of randomly situ- 
ated and randomly oriented strands of ~ 30Athickness. For 98% aerogel an 
estimated average distance between the strands ~ 200Ais of the order of the 
superfluid correlation length £ , which varies in 3 He depending on a pressure 
in the interval 160 — 500A. Estimated mean free path / ~ 1500 — 1800Ais 
greater then £ . According to the theory of superconducting alloys [2] at the 
unconventional Cooper pairing even nonmagnetic impurities lower a tem- 
perature of transition in the superfluid (or superconducting) state T c to the 
extent of a ratio £o/l [3]. Such suppression is observed in 3 He. Below T c in 
a system 3 He+aerogel two superfluid phases are observed [4]. By the anal- 
ogy with the bulk 3 He one of the observed phases is referred as A-like phase 
and the other - as B-like. Recent observation of homogeneously precessing 
domain in the B-like phase [5] is a very strong evidence that the order pa- 
rameter of the B-like phase is close to that of BW-phase of the bulk 3 He. 
Situation with the A-like phase is not so clear. Identification of that phase 
rises important question, whether in a presence of aerogel new superfluid 
phases of 3 He can realize, with the order parameter which is different from 
that of possible superfluid phases of bulk 3 He. Answer to this question is a 
main objective of the present paper. In what follows a procedure for deter- 
mining the order parameters of superfluid phases 3 He in aerogel in a vicinity 
of T c is formulated and applied to equal spin pairing (ESP) states which are 
good candidates for description of the A-like phase. 



2 Interaction of aerogel with the order pa- 
rameter 

At Cooper pairing with / ^ except for the above mentioned general sup- 
pression of T c effect of fluctuations in positions and orientations of the strands 
of aerogel is of importance. In a vicinity of T c this effect can be described 
phenomenologically by a random field, acting on the order parameter. Cor- 
responding term in the free energy can be constructed using symmetry ar- 
gument. In 3 He Copper pairs are formed in a state with the orbital moment 
/ = 1 and spin s = 1. In that case the order parameter is a 3x3 complex 
matrix A^j, index \x refers to spin and j to orbital projection. Interaction of 
3 He with the strands of aerogel is due to the scattering of quasiparticles on 
the strands. Quasiparticles change their momenta and interact directly with 



2 



the orbital part of the order parameter. Aerogel consists of magnetically in- 
ert SiC>2, but when immersed in the liquid 3 He it forms a layer of "solid" 3 He 
on the strands. Quasiparticles of the liquid can exchange their spins with 
that of atoms in the "solid" layer of 3 He. This exchange provides a mech- 
anism of interaction of aerogel with the spin part of the order parameter. 
In experiments this mechanism can be switched off by addition of 4 He in a 
cell. Atoms of 4 He substitute 3 He in the "solid" layer and make exchange 
impossible. Admixture of 4 He strongly influences magnetic properties of su- 
perfluid 3 He in aerogel. In what follows it is assumed that the strands of 
aerogel are covered by a layer of 4 He. In that case aerogel interacts only with 
the orbital part of the order parameter A^j and in a principal order on 
a corresponding contribution to the free energy has the form [6]: 

F„ = N(0) J vMA^A^r, (1) 

where N(0) is the density of states at the Fermi level and 77^ (r) - a random 
static tensor field. On the strength of t — > — t invariance tensor rjji(r) is 
real and symmetrical, its isotropic part ^r)u(r)8ji describes local variations of 
T c = T c (r) due to fluctuations of the density of scatterers. Anisotropic part 
r]ji(r) — ^r]u(r)5ji = r/jf describes local splitting of T c for different projections 
of angular momenta because of the breaking of spherical symmetry by the 
aerogel strands. Isotropic part of the random field will be included in T c = 
T c (r). Absolute value of the random field can be estimated with the use of the 
results of the paper [7]. By the order of magnitude it is ~ (xC,q/R) ~ 
(£o/0) "where / - is a mean free path, R- average radius of a strand x - 
fraction of a volume, occupied by aerogel. For a 98% aerogel £oA ~ 1/10. 
The field rjji(r) varies on a scale of the average distance between the strands 
d ~ Rj ^fx. In a 98% aerogel it is of the order of £o- If the order parameter 
varies on a scale of ~ d a loss of the gradient energy ~ {C,o/d) 2 exceeds a 
gain of energy in the field r]ji(r) to an extent of ^ ^> 1. For that reason 
the order parameter does not follow variations of the field and does not form 
states localized on a scale ~ d. It is possible thou to form localized states 
with a scale L^> d even in a weak field %(r) because of a degeneracy of the 
average order parameter A^j over rotations in the orbital (momentum) space. 
According to Imry and Ma [8] a random field can destroy a long range order 
if the ordering is characterised by a continuously degenerate order parameter. 
They show in particular that for a vector order parameter s(r) interacting 
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with a random field h(r) as: 

F IM = -J s(r) • h(r)d\ (2) 

long range order is destroyed by arbitrary small field. The argument is follow- 
ing. The average value of the random field h(r) is zero, i.e. -p- / h(r)<i 3 r — > 
at L — > oo, where L-is a linear scale of a region of integration. The integral 
tends to zero as ~ (d/L) 3 ^ 2 . As the same power of L tends to zero a gain 
of energy of the order parameter due to its orientation in a field averaged 
over a region with a scale ~ L. The loss of gradient energy decreases faster 
~ (£ /L) 2 and for large L formation of domains becomes advantageous. As 
a result the long range order is destroyed. General argument of Imry and 
Ma does not apply directly to the superfluid 3 He in aerogel. The interaction 
F v differs from Fjm in a sense that there exist finite A^ for which F v is zero 
for all possible 77^. These A^ can be found from the equation 

vfKiKx = 0, (3) 

or, equivalently from the equation which does not contain 77^: 

KiKj + A ri A U = 5 H ' const ( 4 ) 

This equation determines the real part of a product A^A*^, its imaginary 
part is an arbitrary antisymmetric tensor. When interaction with a random 
field turns to zero local reorientation of the order parameter does not give a 
gain in the energy and a long range order is preserved. From that qualita- 
tive argument one concludes that in a situation when average value of the 
order parameter gives a principal contribution to free energy eq. (3) gives a 
necessary condition of stability of a phase corresponding to A^ with respect 
to a random field r/ji(r) [9]. In the next section a procedure of finding of the 
order parameter A^ in a presence of a random field 7?j/(r) is formulated. It 
will be shown that the criterion (3) or (4) emerges as a natural requirement 
of consistency of the procedure of minimization of the corresponding free 
energy. 
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3 Selection of superfluid phases 



With the account of interaction eq.(l) Ginzburg-Landau functional takes the 
following form: 

f gl = mo) [ dV[Tv; + vjii^A^ + fois+ 

J Z 8 = 1 

- ^—^—^ + K 2 —^—^ + K 3 —^—^ , 5 

2 \ oxj oxj oxj oxi oxj ox i y J 

where r = (T — T c )/T c , I s - 4-th order invariants in the expansion of the free 
energy over A^j, we don't need their explicit expressions (cf. [10]). Coeffi- 
cients /3i, ...05, Ki, K2, K3 - phenomenological constants. In what follows we 
use for their evaluation the weak coupling values, in particular we assume 
Ki = K 2 = K 3 = K. The gradient terms can acquire random corrections as 

OA* 

well e.g. of a form: Uj{v)A^i^ L 1 where Uj(r) is a random vector . Up to a 
prefactor Ti/m this is a local random velocity. The terms of such form are ob- 
tained by extension of derivatives ^- — > -£^ + Uj(r) in the expression for free 

energy (5). These terms do not influence the selection of superfluid phases, 
that's why they are not kept in what follows. Variation of the functional (5) 
over A*^ renders an equation for the equilibrium order parameter: 

and variation over A^j - the complex conjugated equation. 

Random field rjji(r) according to the above estimation is small. Effect 
of a small random field on a conventional superconductor in a vicinity of T c 
was analysed by Larkin and Ovchinnikov [1 1] . More complicated form of the 
order parameter in 3 He and its degeneracy requires a nontrivial modification 
of that procedure. 

Random field induces fluctuations of the order parameter a w - at its aver- 
age value Apj, i.e. -A w -(r) = A^ + a w (r). A long range order is established 
when A^j 7^ 0. This condition determines the transition temperature T c . Not 
too close to T c fluctuation a w - is magnitude of the first order on rjji. Let us 

1 V.I.Marchenko pointed out to me a possibility of their existence 
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expand eq. (6) at = keeping terms up to the second order on and 
Vji- 

1 5 dl s d 2 I s d 2 I s 



2 \dA* 4 dA m dA m Pl dAr.MT m dAfn 



1 T , / 9 2 a u1 - 9 2 a u ; 
-if — + 2- 



(7) 

Now let us take average of the obtained equation over the length scales, which 
are much greater then the average distance between the strands of aerogel. 

1 5 dl 1 d 3 I 

TA »' + 2 5 M dt + ^dA U dA md A m < ""'> a » > + 

d 3 I x 

< al n a P i >)] = -< a/dTjij > . (8) 



8A*„dA*dA 



Except for A^ the obtained equation contains averages of binary products 
of fluctuations < a vn api > etc.. Equations for are obtained by separation 
of fast varying terms in eq.(7) and in its complex conjugated 

1 5 d 2 I s d 2 I s 

1 5 (9 2 / (9 2 / 

w 2 s=1 dA^jdA* n dA^jdAyn 

H^ +2 29'^*' (io) 

This is a linear inhomogeneous system of equations. As a consequence of 
the discussed above degeneracy of the order parameter A^j the corresponding 



6 



homogeneous system has nontrivial solutions. These are variations of and 
.A* ■ at an infinitesimal rotation 6 n : 

= O n e jnr A^ r , u; 3 = 9 n e> nr A; r , (11) 

where e j,nr - the absolutely anti-symmetric tensor. It is convenient to use 
Fourier transformed equations (9), (10) for 77^ (k) and a w (k). For the further 
argument only a character of singularity of a w -(k) at k — > is important. 
To find it we neglect anisotropy of the gradient terms in eqns.(9),(10) and 

substitute \K (*2ji) and \K (^) instead of \K + 2^) and 

1 K( d2a l j I 2 ^ ^ 

2- fv dxidxjj- 

Taking product of eq. (9) by a;*-, eq. (10) by cj w - and taking sum of 
both equations one arrives at the following expression for the projection 

a W ( k Kj + «( k V« = °"( k ) : 

^(k) = -^(^g ri r // ( ; ) (k)), (12) 

where Q T \ = A^A*^ + A^A*^ . This projection gives a contribution to the 
averages < a vn api > which is proportional to 

J < a w (k)a w (-k) > d 3 k ~ 

9 P e^Q rl 9 q e m « s Q sn J < ^(k)ifi(-k) > ^. 

Form of the correlation function in the integrand after averaging over direc- 
tion of k is determined by the symmetry properties of r)ff: 

J < i#>(k)i#>(-k) > ^ = S(fc)[2M™ - 3(Vmn + (13) 
Finally 

y < a^(k)a w (-k) > d 3 k ~ 

^ 6 ^ QrlQsn ~t~ (^pt/^rs ^ps^gTO^rn^sn] ^ ^2 (^^O 

The integral in the left-hand side diverges at the lower limit and the averaged 
products < a w a w > give diverging contributions to eq.(5). To make the 
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iteration scheme eqns.(5),(6) consistent one has to require the coefficient in 
front of the diverging integral to be zero. Equating to zero the expression 
in square brackets in eq. (11) one arrives at Q rn = q5 rn , where q is a real 
number. In terms of the order parameter this condition reads as: 

A vr A U j -\- A U jA ur = const ■ 5 r j. 

that coincides with the criterion (4). 

The order parameters meeting criterion (4) can be refered as "quasi 
isotropic" since the energy of their interaction with aerogel does not change 
at arbitrary rotation in the orbital space, or continuous degeneracy is pre- 
served in a presence of the random field i]ji(r). The order parameter A^j 
enters the tensor of superfluid densities in the combination A^A* ^ + A^A*^ 
i.e. condition (4) is the requirement of the isotropy of that tensor. It might 
be well to point out here the difference between the argument of Imry and 
Ma and the argument presented in this section. In the former case the order 
parameter is assumed to be fixed by minimization of the free energy without 
the random field. Possible disordering is caused by a random walk of the 
order parameter over its space of degeneracy. No other change of the order 
parameter is admitted. In the latter case all variations of the order param- 
eter A^j are considered. That makes possible an adjustment of the order 
parameter to the random field. To emphasize the property of phases, which 
satisfy criterion (12) or (13) to resist disordering action of the random field 
a term "robust" was suggested [12]. 

Hence, the procedure of finding of order parameters of possible superfluid 
phases of liquid 3 He in a presence of the random field ^(r) begins with a 
choice of a family of matrices A^ meeting criterion (4). These matrices form 
a "zero-order approximation" for the order parameter. As a next step the 
first order corrections and a*^ are expressed in terms of A^ and ^(r) 
via eqns. (9) and (10). In a practice it is more convenient to work with the 
Fourier transforms a w -(k). After calculation of the averages < a un api > etc. 
and substitution of these averages in eq.(8) it becomes a closed equation for 
A^j. Coefficients Pi,...p^,K and correlation functions < r) un (k.)r)(3i(—'k) > 
have to be given. When rjji(r) = a conventional equation for the extrema 
of free energy of pure 3 He is recovered. 

The order parameter of BW-phase 

A™ = Ae*i2 w -, (15) 
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where R^j - a real orthogonal matrix, satisfies criterion (4). If the dipole 
interaction is neglected this matrix up to a relative rotation of spin space 
with respect to the orbital is equivalent to the unit matrix 8^. Aerogel 
is assumed to be uniform and isotropic. In that case tensor structure of 
correlation functions < 77^ (k)^^— k) > is determined by symmetry [6]. It 
is clear that with that functions the order parameter proportional to the 
unit matrix satisfies eq. (8) and the BW-phase can be a minimum of free 
energy in a presence of aerogel. In comparison with the pure 3 He the values 
of phenomenological coefficients /3i, ..f3§ will change. This influences a region 
of stability of BW-phase and changes thermodynamic quantities, depending 
on these coefficients. We don't consider here a problem of explicit calculation 
of corrections to Pi, ...fl 5 . 

The order parameter of ABM-phase 

A ri = A ^d fM (m j + ihj), (16) 

does not satisfy criterion (4). This brings up the question: what kind of the 
order parameter can describe the properties of A-like phase. 



4 Zero-order approximation for ESP phases 

The measured magnetic susceptibility of A-like phase coincides with that 
of the normal phase [4]. It means that the order parameter of that phase 
does not contain components corresponding to zero projection of the spin of 
Cooper pair on the direction of magnetic field i.e. it belongs to ESP (equal 
spin pairing) type. General form of the ESP-type order parameter is: 

A m = A ^[d^( m j + irij) + e^lj + ipj)}, (17) 

where and e M are mutually orthogonal unit vectors, mj,rij, lj,Pj - arbitrary 
vectors. Direct substitution of the order parameter (17) into eq. (4) shows 
that the criterion is met if the vectors rrij, rij, lj,Pj satisfy the equation 

rrijmi + rijUi + Ijk + pjpi = 5 jt . (18) 

The normalization condition A^jA*^ = A 2 is used. One of the solutions of eq. 
(18): (p=0, m,n,l- orthonormal basis) has been discussed before [9]. To find 
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all solutions consider four four- dimensional vectors M s , N s , L s , P s (s=l,2,3,4) 
which obey the equation 

M r M s + N r N s + L r L s + P r P s = 5 rs , (19) 

Unique (up to a general rotation and reflections) solution of this equation is a 
set of four orthonormal vectors : ■ q^ = 5 ab . Let us take an arbitrary 
unit four-vector v = (z/ 1; z/ 2 , z/3, z/4) and take projections of the vectors q^ 
on the three-dimensional hyperplane orthogonal to u. As a result we obtain 
four three-dimensional vectors: 

m = q^ — v\v, n = q^ — z/ 2 z>, 1 = q^ — v^v, p = q^ — z/ 4 z>, (20) 

Taking a product of m^mi + rijiii + Ijk + PjPi, on an arbitrary vector ai 
perpendicular to u, and using for m,n,l,p expressions (20) one can see that 
these vectors satisfy eq. (18). With the aid of formulae (20) one can find 



other useful properties of the vectors m,n,l,p: 

m 2 + n 2 + l 2 + p 2 = 3 (21) 

m • n = m • 1 = — ^1^3, n • 1 = -v 2 v 3 , ... (22) 

m 2 = 1 - v\,n 2 = 1 - v 2 ,... (23) 



From eqns. (22) follows that [m x n] • [1 x p] = 0, i.e. the normals to the 
planes defined by the pairs of vectors m,n and l,p are mutually orthogonal. 
This property is preserved for arbitrary choice of pairs among the four vectors 
m,n,l,p. Hence eq. (17) with the vectors m,n,l,p, defined according to eq. 
(20) specifies a three-parameter family of quasi isotropic order parameters of 
ESP-type. Substitution of this order parameters in eq. (5) renders zero-order 
energies of the corresponding phases: 

F (o) A4 

= r A 2 + — [ft + 9ft + ft + 5 (ft + ft) - 4(ft + ft) (vm - ^s) 2 ] • (24) 

Parameters v\, 1S2, V3, enter this expression in a combination A = {i>\v^ — 
V2V3). If (ft + ft) = ft 5 < 0, the minimum of free energy corresponds to 
A = i.e. to 

viv A = 1/2V3. (25) 

In a weak coupling approximation both ft and ft are negative and the in- 
equality (ft + ft) < is well satisfied. Condition (25) has simple physical 
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meaning. The order parameters given by eq. (17) are generally speaking 
nonunitary. The corresponding phases can have a finite spin density, which 
is proportional to e^xA^A*^ or, explicitly to (2A 2 /3)[d x e][n • 1 — m • p]. 
With the aid of eq. (22) one can find that the spin density vanishes if eq. 
(25) is satisfied. This condition specifies two-parametric family of nonferro- 
magnetic quasi isotropic phases which possibly includes A-like phase. A pos- 
sible parametrization of this family is: V\ = sin a sin /3, v 2 = sin a cos /3, u 3 = 
cos a sin /3, z/4 = cos a cos (3. The most symmetrical nonferromagnetic phase 
corresponds to the choice a = ir/4, (5 = ir/4. With these values of parameters 
Vl — y 2 — u 3 — U4 — 1/2 and the absolute values of vectors m,n,l,p are all 
equal to y/3/2. The angles between every two of the four vectors m,n,l,p 
are also equal. Such set is formed by vectors connecting the center of a right 
tetrahedron with its vertices. 

If /3i5 > then A 2 in equilibrium reaches its maximum value 1/4 which 
is reached at v 1 = z/ 4 ; u 2 = —v 3 or at v 1 = — z/ 4 ; v 2 = v 3 . In both cases it is 
the one-parametric family. In the first case it can be parmetrised as v\ = 
1/4 = -^sin7 v 2 = —v 3 = -^cos'j. The most symmetrical ferromagnetic 
phase corresponds to v\ = —1/2, v 2 = v 3 = = 1/2 i.e. it is obtained from 
the most symmetrical nonferromagnetic solution by inversion of one of the 
vectors rrij , rij , I j , pj . 

5 Effect of magnetic field 

In a magnetic field two terms have to be added to free energy. One is 
quadratic in the field: 

f ( H ] = -\x»vH,H u . (26) 

For all ESP phases one of the principal values of a tensor of magnetic suscep- 
tibility Xiiv coincides with the susceptibility of the normal phase Xn- In 
a vicinity of T c tensor Xfiv is determined by symmetry argument x^v = 
XnSfiv — K-iAnjAlj + A U jA*j). Second term in the r.h.s. describes a decrease 
of the transverse susceptibility in comparison with the normal value Xn- This 
is two-dimensional tensor with the principal values ^-Ai^, where Ai^ are 
the roots of the equation 

A 2 - 3A + 2 + A 2 = 0. 
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In a nonferromagnetic phase A = 0, Ai = 2, A 2 = 1 i.e. transverse suscep- 
tibility is anisotropic. In a ferromagnetic phase A 2 = 1/4, Ai = A 2 = 1, 
and transverse susceptibility is isotropic. Equilibrium orientation of the or- 
der parameter corresponds to maximum susceptibility Xnv in the direction of 
magnetic field and additional energy eq. (26) is the same for all ESP-phases. 

Except for the quadratic there is a linear on magnetic field term in the 
free energy: 

ff = iCe^xA^AljHx. (27) 

In the bulk 3 He this term gives rise to the splitting of the transition from 
the normal phase in the A-phase in two transition with close temperatures. 
Ferromagnetic A-phase appears first, it contains Cooper pairs with only 
one spin projection. At a lower temperature transition in A 2 -phase occurs 
with both spin projections present in the condensate. The coefficient ( is 
proportional to the derivative of the density of states over energy and the 
temperature interval where Ai-phase exists is small as /j,H/ef, where fi is 
the magnetic moment of 3 He nucleus and Bp - Fermi energy. 

Let us take into account additional free energy eq. (27) in a presence of 
aerogel. 

p(") /- TT A o A 4 A 4 

^gy = (r ~ ^)A 2 - — A 5 A 2 + — [ft + 9ft + ft + 5(ft + ft)]. (28) 

This expression has to be minimized over A and over A 2 . The result of 
minimization depends on a sign of ft.5. If fts > 0, then for all A 2 the 
energy reaches its minimum at |A| = 1/2, i.e. the ferromagnetic phase is 
stable. Transition in the superfluid state takes place at r = (H/6. At 
r < (H/6 A 2 = -9r H /B, where r H = r - (H/6, B = 9ft + ft + 5ft + 4ft. 
Magnetic moment has a field independent increment which is proportional 
to A 2 : M = N(0)CA 2 /6. 

If fts < 0, the ferromagnetic phase minimizes free energy (28) only 
in the temperature interval [1 + (B / p 15 )]((H/6) < r < (H/6. At r 2 = 
[1 + (B/Pi 5 )]((H/6) a second transition takes place in the phase with A = 
— 4 p^ A 2 ■ When temperature is lowered A — ► i.e. an additional magnetic 
moment vanishes. The transition at r = r 2 is analogous to A\ — > A 2 transi- 
tion in bulk 3 He. Hence the considered ferromagnetic phases are analogous 
to Ax-phase of bulk 3 He except that for the new phases the pairing amplitude 
is finite for both projections of spin s = 1 and s = — 1. 
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6 Discussion 



We conclude that in 3 He in a presence of aerogel only such superfluid phases 
can exist for which the order parameter satisfies criterion (4). The examples 
of A-like and Ai-like phases demonstrate that this condition does not deter- 
mine a matrix uniquely, but selects a family of such matrices. For the 
further selection a next approximation on the random field 77^7 (r) has to be 
considered. The answer depends on correlation functions of the random field, 
which are not known. A class of admitted solutions can be restricted with 
the aid of the observed properties of the phases. For example the splitting 
of the phase transition in a magnetic field is the evidence that /3i5 < and 
away from the T c nonerromagnetic phase is stable. If the A\ — > A 2 transi- 
tion is absent then a ferromagnetic phase is stable (/3i5 > 0). Experimental 
situation in that respect is not quite clear. 

A characteristic property of all quasi isotropic phases, which can be used 
for a general check of the proposed scheme is the isotropy in a principal order 
on r)ji(r) of the tensor of superfluid densities. In particular this is a difference 
between the proposed A-like phases and of the A-phase of bulk 3 He. 
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